Theory of Model Kohn-Sham Potentials and its Applications

The purpose of Kohn-Sham density functional theory is to develop increasingly accurate approximations to the exchange-correlation functional or to the corresponding potential. When one chooses to approximate the potential, the resulting model must be integrable, that is, a functional derivative of some density functional. Non-integrable potentials produce unphysical results such as energies that are not translationally or rotationally invariant. The thesis introduces methods for constructing integrable model potentials, developing properly invariant energy functionals from model potentials, and designing model potentials that yield accurate electronic excitation energies. Integrable potentials can be constructed using powerful analytic integrability conditions derived in this work. Alternatively, integrable potentials can be developed using the knowledge about the analytic structure of functional derivatives. When these two approaches are applied to the model potential of van Leeuwen and Baerends (which is non-integrable), they produce an exchange potential that has a parent functional and yields accurate energies. It is also shown that model potentials can be used to develop new energy functionals by the line-integration technique. When a model potential is not a functional derivative, the line integral depends on the choice of the integration path. By integrating the model potential of van Leeuwen and Baerends along the path of magnitude-scaled density, an accurate and properly invariant exchange functional is developed. Finally, a simple method to improve exchange-correlation potentials obtained from standard density-functional approximations is proposed. This method is based on the observation that an approximate Kohn-Sham potential of a fractionally ionized system is a better representation of the exact potential than the approximate Kohn-Sham potential of the corresponding neutral system. Removing 1/2 of an electron leads to the greatest improvement of the highest occupied molecular orbital energy, which explains why the Slater transition state method works well for predicting ionization energies. Removing about 1/4 of an electron improves orbital energy gaps and, when used in time-dependent density functional calculations, reduces errors of Rydberg excitation energies by almost an order of magnitude

Electron affinity of liquid water.

Understanding redox and photochemical reactions in aqueous environments requires a precise knowledge of the ionization potential and electron affinity of liquid water. The former has been measured, but not the latter. We predict the electron affinity of liquid water and of its surface from first principles, coupling path-integral molecular dynamics with ab initio potentials, and many-body perturbation theory. Our results for the surface (0.8 eV) agree well with recent pump-probe spectroscopy measurements on amorphous ice. Those for the bulk (0.1-0.3 eV) differ from several estimates adopted in the literature, which we critically revisit. We show that the ionization potential of the bulk and surface are almost identical; instead their electron affinities differ substantially, with the conduction band edge of the surface much deeper in energy than that of the bulk. We also discuss the significant impact of nuclear quantum effects on the fundamental gap and band edges of the liquid

Local and Global Effects of Dissolved Sodium Chloride on the Structure of Water

Determining how the structure of water is modified by the presence of salts is instrumental to understanding the solvation of biomolecules and, in general, the role played by salts in biochemical processes. However, the extent of hydrogen bonding disruption induced by salts remains controversial. We performed extensive first-principles simulations of solutions of a simple salt (NaCl) and found that, while the cation does not significantly change the structure of water beyond the first solvation shell, the anion has a further reaching effect, modifying the hydrogen-bond network even outside its second solvation shell. We found that a distinctive fingerprint of hydrogen bonding modification is the change in polarizability of water molecules. Molecular dipole moments are instead insensitive probes of long-range modifications induced by Na⁺ and Cl^– ions. Though noticeable, the long-range effect of Cl^– is expected to be too weak to affect solubility of large biomolecules

Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals

We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice

Density and Compressibility of Liquid Water and Ice from First-Principles Simulations with Hybrid Functionals

We determined the equilibrium density and compressibility of water and ice from first-principles molecular dynamics simulations using gradient-corrected (PBE) and hybrid (PBE0) functionals. Both functionals predicted the density of ice to be larger than that of water, by 15 (PBE) and 35% (PBE0). The PBE0 functional yielded a lower density of both ice and water with respect to PBE, leading to better agreement with experiment for ice but not for liquid water. Approximate inclusion of dispersion interactions on computed molecular-dynamics trajectories led to a substantial improvement of the PBE0 results for the density of liquid water, which, however, resulted to be slightly lower than that of ice

Photoelectron Spectra of Aqueous Solutions from First Principles

: We present a combined computational and experimental study of the photoelectron spectrum of a simple aqueous solution of NaCl. Measurements were conducted on microjets, and first-principles calculations were performed using hybrid functionals and many-body perturbation theory at the G0W0 level, starting with wave functions computed in ab initio molecular dynamics simulations. We show excellent agreement between theory and experiments for the positions of both the solute and solvent excitation energies on an absolute energy scale and for peak intensities. The best comparison was obtained using wave functions obtained with dielectric-dependent self-consistent and range-separated hybrid functionals. Our computational protocol opens the way to accurate, predictive calculations of the electronic properties of electrolytes, of interest to a variety of energy problems

First-Principles Simulations of Liquid Water Using a Dielectric-Dependent Hybrid Functional

We carried out first-principles simulations of liquid water under ambient conditions using a dielectric-dependent hybrid functional, where the fraction of exact exchange is set equal to the inverse of the high-frequency dielectric constant of the liquid. We found excellent agreement with experiment for the oxygen–oxygen partial correlation function at the experimental equilibrium density and 311 ± 3 K. Other structural and dynamical properties, such as the diffusion coefficient, molecular dipole moments, and vibrational spectra, are also in good agreement with experiment. Our results, together with previous findings on electronic properties of the liquid with the same functional, show that the dielectric-dependent hybrid functional accurately describes both the structural and electronic properties of liquid water